The Optimal Faucet 
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The production of small fluid droplets relies on an instability of solutions to the Young-Laplace 
equation. We ask whether smaller droplets can be produced by changing the shape of the nozzle. 
At a given critical pressure, the circular nozzle actually produces the largest droplet. The droplet 
volume can be decreased by up to 18% using a triangular nozzle with stretched corners. 

PACS numbers: 68.03.Cd, 47.20.Dr, 02.30.Xx 



A standard protocol for producing small droplets is as 
follows: a pipette, of circular cross-section, is pressurized 
at one end, pushing out a small fluid droplet. If the noz- 
zle is sufficiently small, force balance requires that the 
droplet has constant mean curvature. At a critical pres- 
sure, this equilibrium shape becomes unstable, ultimately 
leading to the droplet detaching from the nozzle. 

The volume of fluid entrained during this process is 
set by the total fluid volume contained in the critical 
droplet. This volume scales like r 3 , where r is the noz- 
zle radius. On the other hand, the critical pressure for 
ejecting this droplet scales like 7/r, where 7 is the liquid 
surface tension. Thus, ejecting smaller droplets requires 
higher pressures. The smallest size droplet that can be 
ejected is thus determined by the highest pressure that 
can be reliably applied to the nozzle, without material 
failure, etc. 

One strategy for creating smaller droplets than those 
dictated by the instability of a static droplet is to use 
a time varying forcing at the nozzle. This method has 
achieved an order of magnitude decrease in droplet vol- 
ume 0- 

However, typical nozzles use a circular cross section. It 
is not unreasonable to imagine that changing the shape 
of the cross section to be some other shape may decrease 
the ejected droplet volume, while maintaining the same 
applied pressure. For example, imagine that we have 
a circular nozzle with a pendant droplet just below the 
critical volume: by "squeezing" the shape of the nozzle 
cross section into an elliptical shape, one might cause the 
droplet to detach at a lower volume. 

In this paper we address the question: what is the 
shape of a nozzle for which the ejected droplet volume is 
minimized, for a given applied pressure? We demonstrate 
that circular nozzles do not eject the smallest droplets; in- 
stead, the optimal nozzle more closely resembles an equi- 
lateral triangle, albeit with "stretched" corners. The best 
nozzle shape that we have found has an ejected droplet 
volume about twenty percent smaller than the circular 
nozzle with the same critical pressure. Our method is 
inspired by and extends J. Keller's classic treatment of 
the Euler buckling problem with a beam of nonuniform 



cross section 0. Recently, the method has been applied 
to the optimization of a bistable switch |3| • For a detailed 
mathematical treatment of capillary surfaces in general, 
see 0. 

This Letter is organized as follows. We first explain 
the origin of the pendant droplet instability. Then we 
describe our method for reducing droplet size. Lastly we 
provide numerical calculations implementing the method, 
and present the candidate optimal nozzle. 

Pendant Droplet Instability. - The instability of a 
droplet protruding from a nozzle is due to a bifurcation, 
most easily seen in the case of a circular nozzle that is 
much smaller than the capillary length, which allows us 
to neglect gravity. The shape of the droplet is then de- 
termined by the Young-Laplace equation p = jK, where 
p is the pressure difference across the liquid/air interface, 
7 is the surface tension, and K is the mean curvature of 
the droplet surface. This equation describes a surface of 
constant mean curvature p/j with the nozzle edge as its 
boundary. If the boundary is a circle, then the solution 
must be a section of the sphere with mean curvature p/j. 
From the familiar relation 
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(1) 



we deduce that the radius of curvature of the droplet is 
2'y/p. For small p, such that the sphere radius is much 
greater than the nozzle radius, the solution is a shallow 
spherical cap. But note that its complement, the rest of 
the sphere, is also a solution. As p is increased, these 
two solutions approach each other until both become a 
hemisphere with the nozzle at the equator. The pressure 
at which the two solutions meet is the critical pressure 
p* , and the corresponding degenerate solution is unsta- 
ble. Note that the critical pressure is also the maximum 
pressure, for the nozzle cannot support a sphere smaller 
than itself. 

For a noncircular nozzle, we no longer have such a sim- 
ple geometric picture, however key features remain. The 
unstable solution is still characterized by a bifurcation 
at which two solutions meet, corresponding to the maxi- 
mum pressure achievable for the given nozzle. The criti- 
cal pressure for a general nozzle can be computed as fol- 
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lows: let the droplet surface be parameterized as a func- 
tion R(u, v) over a domain D in the uu-plane, which takes 
value in three dimensional physical space. The bound- 
ary of the domain dD corresponds to a closed curve C 
which represents the nozzle. The curvature is a nonlinear 
functional of the surface and its derivatives up to second 
order, hence the equation for the droplet shape has the 
form 

jK[R,VR,VVR]=p, (2) 

where V is the gradient operator in the uv-pl&ne. 

Upon increasing the pressure p — ► p + Sp, the surface 
changes: R — > R + 5R. Equation J5J implies that the 
variation 6R and Sp are related by 



jLSR = Sp, 



(3) 



where LSR is the change in mean curvature induced by 
the surface change. L is a differential operator acting on 
SR. 

At the critical solution, the pressure is at a maximum; 
therefore, there must be a solution w = SR to equation 
® with Sp = 0. The solution w satisfies 



Lw = 



(4) 



with boundary condition w = at dD. Note that the 
pressure dependence in this formula arises because L = 
L[R] depends implicitly on the pressure p through R. 
Hence, the existence of a nonzero it; is a diagnostic for 
finding the critical solution to J2J and the corresponding 
critical pressure p* . 

Optimization Method. - Now, to find the optimal noz- 
zle, we need to derive a relation between the change in 
critical pressure and change in nozzle shape. Since pres- 
sure and volume are conjugate variables, increasing criti- 
cal pressure is tantamount to decreasing critical volume. 
By iteratively changing the nozzle shape to increase crit- 
ical pressure, we will thus arrive at a nozzle which pro- 
duces smaller droplets. We compare the critical volume 
of the deformed nozzle with that of the circular nozzle 
that corresponds to the same critical pressure, since pres- 
sure is the control variable in practical situations. 

Suppose that a given nozzle shape C has a critical pres- 
sure p* , a critical droplet shape R* , and a corresponding 
w. All of these quantities change when the nozzle shape 
C — > C + SC. The change in the droplet shape SR is 
linearly related to the pressure change Sp by equation 
© with the boundary condition SR = SC at dD. On 
the other hand, since the critical solution maximizes the 
critical pressure, w does not change to leading order in 
SC. 

The change in critical pressure induced by SC can 
therefore be computed by taking the inner product of 
both sides of (J3J with w: 

l(w,L6R) = j(5R,Lw) + 7 / b{SR,w) 



= + 7^) b(SC,w) 
= (w,Sp). 



Therefore 



Sp = 



1 §b{SC,w) 

M 



(5) 



Here 6(», •) denotes the boundary integrand from inte- 
grating by parts. The derivation also uses the self ad- 
joint ncss of L, which is readily demonstrable by explicit 
computation y|. Equation Q is an explicit relation be- 
tween a change in the nozzle shape (SC) and the resulting 
change in critical pressure. 

Explicit Formula for Sp. - We choose the nozzle C 
to lie in the xy-plane, enclosing the origin. Then the 
droplet surface may be given in spherical coordinates by 
the distance from the origin (R) as a function of the two 
angles 9 6 [0,7r/2] and <f> E (0, 27r]. To avoid the coor- 
dinate singularity at the pole (9 — 0) we use u, v given 
by u = tan(0/2) cos(0) and v = tan(0/2) sin(<^>). Hence 
the surface is a scalar function R(u,v); its domain D is 
the unit disk in the uw-plane. We retain <fi to denote the 
polar angle in the uu-plane. 

An appealing feature of this coordinate system is that 
the line element remains diagonal: 



ds 2 



dR 2 + T(du 2 + dv 2 ), 



where T — 4R 2 /(1 + u 2 + v 2 ) 2 . It is then straightforward 
to compute the free energy E = J (jdA — pdV) which 
yields, upon variation, the Young-Laplace equation 



V • (CVR) + AR = F, 



(6) 



where V is the usual gradient operator in the uw-plane. 
The coefficients are 



C = 



1 



A = C 



F = p 




(l+P 2 ) 



plane. dD corresponds to p = 1. 

In our coordinate system, the pressure change is 



Sp 



S W V 



d<t>8C 



w p R{R 2 + R 2 ) 
(R 2 +R 2 „ + Rlf/ 2 ' 



(7) 



where S W V = J d 2 p w n^?p ■ Here and in the following 
we use subscripts to denote partial differentiation. 

We can recast this expression into a form that is 
more geometric. First, the contact angle a between the 
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drop and the plane of the nozzle is given by cota(0) = 
R P /(R 2 + R 2 ( f > ) 1 l 2 \ p =i where the right hand side is evalu- 
ated at the boundary. Second, we define w± = w p / (R 2 + 
R^ > ) 1 l 2 \p=\ which can be understood as follows - note 
that w is the difference between the outer and inner so- 
lutions as the pressure approaches bifurcation. Using the 
contact angle given above, this expression is the differ- 
ence between the slopes (with respect to the vertical) of 
the outer and inner solutions at the boundary. This is a 
coordinate independent quantity. Third, we observe that 



d<j>8cR = 




= dl SN, 



where dl is the line element, and SN is the change of 
the nozzle in the direction locally normal to the nozzle. 
Lastly, the denominator 5 W V in J7J) is just the change in 
volume from changing the surface by w. Putting these 
facts together, the pressure change is 

Sp = \- <L dl SN w± sin 3 a, 
OwV J 

which leads to the prescription for changing the nozzle 

SN ~ -J— w± sin 3 a. (8) 
S W V 

Clearly, for the circular nozzle, symmetry implies that 
SN should be constant. But this amounts to a mere 
reduction in the size of the nozzle; the shape remains a 
circle. So the circular nozzle is at an extremum, in fact 
a minimum of critical pressure for fixed nozzle area. 

For a noncircular nozzle, the contact angle isn't con- 
stant, and hence the change according to the above for- 
mula cannot be constant. So one may change the critical 
pressure while fixing the nozzle area. Moreover, since the 
circular nozzle is the only one (except the infinite strip) 
with a constant contact angle, the process of deformation 
does not end. 

We apply © iteratively to a perturbed circular noz- 
zle to see how the shape evolves away from the circle. 
Figure ^ shows the result of iterations starting with a 
circle deformed by a perturbation with a three-fold sym- 
metry. The perturbation grows with each iteration, and 
eventually the nozzle shape becomes concave. With each 
iteration, we have applied a rescaling in order to maintain 
the nozzle area. Without the area constraint, the nozzle 
would become arbitrarily small in accordance with J^J. 
We are interested in the shape of the nozzle, not its size. 
We also apply the Savitzky-Golay filter [j| at each itera- 
tion to smooth out the mesh noise. The solutions to the 
Young-Laplace equations are obtained using the nonlin- 
ear PDE solver in the MATLAB® PDE Toolbox, which 
implements the finite element method for elliptic equa- 
tions with variable coefficients, exactly of the form in 10. 
For each nozzle shape, we start at a pressure below the 




FIG. 1: Evolution of nozzle shape with threefold symmetry, 
a) Initial nozzle: V = 1.00; b) V = 0.97; c) V = 0.88; d) 
V = 0.82. V is a normalized volume given by |0. 



bifurcation and by choosing different trial solutions ob- 
tain both solutions. Then we bring both solutions to just 
below the critical pressure by stepping up the pressure, 
using the solution at each step as the trial solution for the 
next step. We then use the average of the two solutions 
for our surface, and their difference for w. The validity 
of this procedure can be rigorously shown for a circular 
nozzle, and we expect it to remaind valid for noncircular 
nozzles as long as the pressure is brought close to critical. 

In order to compare and select among nozzle shapes, 
we need a measure of optimality independent of size. 
For every nozzle, we rescale its critical volume by the 
critical volume corresponding to the circular nozzle with 
the same critical pressure. This dimensionless volume is 
given by 




Figure [21 shows a particular sequence of critical prop- 
erties obtained through our iteration procedure. We 
see that the critical pressure begins to increase rapidly 
about the fifth iteration, after which the decrease in V 
slows down, and the nozzle shape becomes stretched out 
(see Figure ^| . This means that in order to decrease 
droplet size at a given pressure, one should use a nozzle 
shape that is roughly triangular, perhaps with somewhat 
stretched out corners; but further deformation does not 
lead to significant improvement. Moreover, gravitational 
instabilities will inevitably become relevant if the "arms" 
become too long 0, • 

It should be emphasized that we have shown a particu- 
lar example of an improved nozzle, generated by a choice 
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FIG. 2: Sequence of iterations away from the circular nozzle 
with an initial three-fold perturbation. The normalized crit- 
ical volume given by © is shown in the bottom graph. The 
arrows indicate the corresponding shapes in Figure Q 



of the initial perturbation. We have tried other pertur- 
bations, leading to shapes with, say, four-fold symmetry 
or without any symmetry, but the three-fold perturba- 
tion has yielded the biggest reduction in the normalized 
critical volume. 

So far we have ignored the effects of gravity, but our 
formalism applies just as well to the problem with grav- 



ity. Including gravity means that the pressure would 
no longer be constant throughout the drop surface, but 
rather a linear function of height: p — > p — p m gh(u,v), 
where p m is the mass density of the liquid, g is the grav- 
itational acceleration, h is the distance below the nozzle, 
and p now denotes the pressure at the nozzle (h = 0). 
Although JBJ acquires a new term as a result, this term 
does not contain derivatives and thus does not contribute 
to the boundary integral. So our formula for the pressure 
change remains the same in the presence of gravity. To be 
sure, the nozzle evolution would differ because the con- 
tact angle and w± will be affected by gravity. Moreover, 
if the nozzle is too large relative to the capillary length, 
then gravity destabilizes all solutions: it is not possible 
to suspend a water drop from a meter wide faucet. It 
would be interesting to examine the case of the interme- 
diate sized nozzle, small enough to have stable solutions, 
yet large enough to be affected by gravity. 
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In computing an inner product (f,Lg) = J dudvfLg, 
where / and g are arbitrary functions of u,v, and the in- 
tegration is over D, we can undo the differentiation on g, 
and instead let the adjoint operator V act on /. Self- 
adjointness (L = U) can be verified explicitly. In the 
derivation, this allows us to simply interchange / and g, 
while introducing the boundary terms from integration by 
parts. 



